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Abstract 

We study CP violation in a multi-Higgs doublet model based on a 53 x Z3 
horizontal symmetry. We consider two mechanisms for CP violation in this 
model: a) CP violation due to complex Yukawa couplings; and b) CP violation 
due to scalar-pseudoscalar mixings. We find that the predictions for e'/e, 
CP violation in B decays and the electric dipole moments of neutron and 
electron are different between these two mechanisms. These predictions are 
also dramatically different from the minimal Standard Model predictions. 
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I. INTRODUCTION 



The origin of CP violation, is one of the outstanding problems of particle physics today. 
So far CP violation has only been observed in the neutral Kaon system. The observed 
CP violation can be explained in many models. It is therefore important to study other CP 
violating experimental observables and compare the results with different model predictions. 
Such study may reveal the real origin of CP violation. 

In the minimal SU(3) C x SU(2) L x U(1) Y Standard Model (MSM), there is only one 
Higgs doublet. When the Higgs develops valcuum expectation value (VEV) v, all fermion 
receive masses. In the mass eigenstate basis, the Higgs coupling to fermions is diagonal, it 
does not mediate CP violating interaction. However, the coupling of the charged current to 
the W-boson becomes complex. It is given by 

L c = ^u i V ij r^ ll d j W+ + H.C. , (1) 

where the matrix Vij is the CKM matrix Vkm [0]- For three generations of quarks, there is 
a non-removable phase in the matrix. This is the source of CP violation in the MSM. This 
matrix is conveniently parametrized as, following Wolfenstein [TJJ 

1 l-A 2 /2 A AX 3 (p-ir])^ 



Vj 



KM 



-A 1 - A 2 /2 - iA 2 X 4 r] AX Z 

i2 



(2) 



I AX 3 (1 - p - irj) -AX 2 1 J 

where A = V us = 0.221. If 77 7^ 0, CP is violated. Unitarity constraints on the matrix 
elements provide very powerful and interesting realtions. The most interesting one is the 
triangle defined by 

v ud v; d + v us v* + v ub v t ; = . (3) 

In the Wolfenstein parametrization, V u d ~ V t b ~ 1 and V ts ~ V* b , we have 

V ub + V t * d « V ub V cb . (4) 
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This defines the triangle shown in Figure 1 with three angles a, (3 and 7. The area of the 
triangle is given by A 2 X 6 r]/2. CP violation in the neutral Kaon system can be explained by 
the "box" interaction ||. If CP violation due to the phase in the CKM matrix is the only 
source for CP violation, experiments at B factories will be able to determine all the three 
angles ||. 

Another class of model for CP violation is the multi-Higgs doublet model. If there are 
more than one Higgs doublets, the neutral scalar couplings to the quarks are not necessarily 
diagonal, and therefore provide new sources for CP violation ||. CP violation can arise in 
three places in this type of models: 1) Non-trivial phase in the Vkm matrix; 2)Non-trivial 
phases in the Yukawa couplings; and 3) Mixings of scalar and pseudoscalar Higgs bosons. 
In cases 2) and 3), CP violation can occurs at the tree level by exchanging neutral Higgs 
bosons. In this talk we will present some studies of CP violation in multi-Higgs models with 
flavour changing neutral currents at the trre level which has CP violation predominantly 
throug mechanisms 2) and 3). 

A most general study suffers from too many free parameters. To have a definite idea, 
we carry out the study in a S3 x Z 3 model proposed by Ma J7J. This model has some 
very interesting predictions for fermion masses and their mixings. It also has interesting 
predictions for CP violation |pHlO|l. We study the predictions for: (i) e'/e; (ii) CP violation 
in the neutral B system; and (iii) the neutron and electron electric dipole moments (EDM). 
We compare these predictions with those in the MSM. 

II. YUKAWA COUPLINGS IN THE S 3 x Z 3 MODEL 

In the S3 x Z 3 model, there are four Higgs doublets, 0i, 2,3,4- The quarks and Higgs bosons 
transform under the S3 x Z3 symmetry as [^] 

<?3L , t R , b R , 0i : (1, 1) , (q 1L , q 2L ) , (03, <Pi) ■ (2, u) , 
{cr, u r ) , (s R , d R ) : (2, u 2 ) , (p 2 ■ (1, w 2 ) , (5) 
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where uj ^ 1, u 3 = 1 is the Z 3 element. The Yukawa couplings consistent with the S3 x Z 3 
symmetry are given by 

L Y = -/l(?lL03«R + q2L<PiC R ) - f 2 q3L<i>ltR ~ hiQlL^R + q2L<p2d R ) 

- U(qxL<fafc R + q 2 L(pAb R ) - f 5 {q 3 L<f>3d R + q^^n) - MsLfabR + E.G. (6) 



where fa = (<$*, —<f)~) T . Without loss of generality we work in a basis where all VEVs 
are real. The up-quark mass matrix is diagonal : M u = Diag{f\v%, /if 4, ftiOi), and the 
down-quark mass matrix can be written as, with a suitable choice of quark phases, 

/ a £b\ 



M° 



(7) 



a b 
V£c c dj 

with a, b, c, d real and £ = |£|e icr complex. M d can be diagonalized by a bi-unitary 
transformation M d = Vi,M d V R . Here M d is the diagonalized down quark mass matrix. Vl 
and V R are unitary matrices. Because the up quark mass matrix is already diagonalized, Vl 
is the CKM matrix Vkm- 

It is convenient to work in a basis of the Higgs bosons in which the Goldstone bosons 
are removed. To this end we define the following |J 



?3 



'•2 



112. 

v 



V2_ 
_ Vl 



V1V4 



V2V4 



1 V4 
\ V 
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— ^ 

1>124 



• -- \ / g \ 



VV\24 

V2V3 
VV124 

V124 
V 

V3V4 

VV124 



H 2 



(8) 



where v\ 2 = v I + v f, v\ 2 a = v i + v \ + u 4j an d f 2 = f 2 + v\ + i>§ + v\. The transformation 
is the same for both the neutral and charged Higgs bosons. For the neutral Higgs bosons, 
G = h° + iGz, where Gz is the Goldstone boson 'eaten' by Z, and h° is a physical field 
whose couplings are the same as the Higgs boson in the MSM. For the charged Higgs bosons 
G is the Goldstone boson 'eaten' by W. In this basis, we have 

L Y = -{D L M d D R + U L M U U R ){\ + —=■) 

vy2 



U L V KM YfD R hf + D L vi M Y?U R hr + E.G. 



(9) 



where hi are the component fields of Hi with Hi = (hf , h® / y/2) . Ul,r = {u, c, t)^ R , and 
Dl,r = {d, s, b)^ R . The couplings % can be easily expressed in term of quark masses, Vl,r, 
and VEVs. 

In general h°' + are not the mass eigenstates. We can parametrize the mixings as 

/ h° \ ( a 00 «oi P'oj \ ( Ro \ 

OikO a ki fi'kj 



Rehl 



\ImhlJ 



\IjJ 



(10) 



where Ri, Ii and i]i are the mass eigenstates, the matrix {a(3) is a 7 x 7 othogonal matrix, 
and (7) is a 3 x 3 unitary matrix. 

The specific values for the mixings depend on the details of the Higgs potential. Unfor- 
tunately they are not determined. To simplify the problem, we will discuss two ) CP 
violation only comes from complex Yukawa couplings; and b) CP violation only comes from 
the mixings of real and imaginary ||10|| . Case a) can be realised by constraining certain 
soft symmetry breaking terms in the potential 0. We further assume, for simplicity, that 
Reh® are the mass eigenstate Ri and consider their effects. The same analysis can be easily 
carried out for I mil® in the same way. The source for CP violation is the non-zero value 
for a which is a free parameter. We will present our results for a = 80°, which is close to 
the maximum of the allowed phase. Case b) can be realised by requiring spontaneous CP 
violation. The value of o will be zero and CP violation arises due to scalar-pesudoscalar 
Higgs boson mixing. For illustration, we consider the effects of a neutral mixed state 



R = cosORehn + sin6Imh\ 
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(11) 



and for the charged Higgs boson we consider mixing 



V + = 722/4 + 723^3 



(12) 



where 7^ are complex numbers, and 1 722 1 2 + 1 723 1 2 = 1- 

The parameters a, b, c, and d are constrained from the down quark masses and the 
CKM mixings. We take as input parameters a = 0.04GeV, b = 0.25GeV, c = 2.66GeV , 
d = 4GeV. The mass eigenvalues for the down quarks are quite insensitive to the phase a. 
For both cases, we have mj, = 4.8GeV , m s = 149MeV and = 9.5MeV. These values are 



well within the allowed regions ||TTJ. The CKM matrix for case a) is 

/ 0.975 -0.222 0.00476 \ 



Vj 



KM 
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0.221 + ^0.0033 0.974 + i0.014 0.043 -^0.0015 

V -0.014 + il.2 x 10- 5 -0.041 -i6.8 x 10- 4 0.998 - i0.034 J 
and for case b) 

/ 0.975 0.22 0.0048 \ 

-0.219 0.975 -0.0436 
V -0.014 0.0415 0.999 J 
The values for the VEV's are not fixed, we only know V3/V4 = m u /m c . We will use 
the values: v\ = v 2 — 4:4GeV, v 3 = 0.9 GeV and t> 4 = 238GeV A for illustration. We shall 
comment on effects of changing these values later. 



KM 
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III. CONSTRAINTS ON THE HIGGS BOSON MASSES FROM THE NEUTRAL 

K AND B MESON SYSTEMS 

The S3 x Z 3 model has very restrictive allowed values for the non-trivial CP violating 
phase in the CKM matrix. The CP violating measure J is less than 2.5 x 10 -6 which is 
too small to explain CP violation in the neutral Kaon sysytem. Therefore in this model CP 
violation due to Higgs boson exchange has to be considered. 

The CP violating parameter e is given by 

ImM* ... 



V2Am K 
6 



where Mf 2 is the matrix element which mixes K° with K°, and Amx is the mass differ- 
ence between m^ h and m^ s . Experimental value for e is 2.3 x lO^e" 1 "/ 4 . The AS* = 2 
Hamiltonian, responsible for M^, generated by exchanging neutral Higgs bosons Ri is given 
by 



We obtain 

jlfg = < K°\H eff \K° > 

;[(«« + ^fcJ^fc 12) + («iM - ^Unii 



fk m K , 5 , ~ d 2 . ~ d ^ 2 



2M|. v 24 (m s + m d ) s 

1 I 2 

+ («fci + ^L)n d i 2 («^ - ia'vJYgni— + - ^ j)2 )) • (17) 

Here we have used the vaccum saturation and factorization approximation results for the 
matrix elements HI3"| . The contribution to the mass difference Amx is given by 2i?eM 12 . 
Similar formula holds for the neutral B system. 

To constrain the Higgs boson masses, we require that the neutral Higgs boson contribu- 
tions to the mass differences in the neutral K and B systems to be less than the experimental 
values: Amx/mx = 7 x 10~ 15 , and Am^/m^ = 8 x 10~ 14 . We find that for case a) the 
tightest constraints on the masses of Reh\ 2 are from the mass difference AM# of the neu- 
tral B mesons which gives > 2.9TeV and Mh 2 > 3.1TeV. With these masses, Reh® 2 
can not produce large enough e. Similar consideration yields M h3 > 3.5TeV, and we find 
the experimental value of e can now be produced if the mass is about 5.6TeV. The mass 
difference AMk of the neutral K mesons gives weaker bounds in all cases. For case b), the 
experimental value of AMg constrains Mr > 3TeV. From the experimental value of e, we 
obtain sinOcosO / 'Mj, = 1.1 x 1(T 8 GeV~ 2 which implies M R < 7TeV. 

IV. PREDICTIONS FOR e'/e. 



In this section we study the direct CP violation in K^ s — » 27r decays. CP violation in 
these processes is characterized by the value of e'/e. e'/e is defined as 
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e^ _ u£ - ImA 2 /ReA 

e ~ £ + ImM 12 /AM K ' [ ' 

where u = ReA 2 /ReA = 1/20, £ = ImA /ReA . Here A and A 2 are the A/ = 1/2, 3/2 
decay amplitudes for K^s — ► 27r. 

In the MSM, the contribution to e'/e is donminantly from the gluon penguin. However, 
for large top quark mass of order 200 GeV, the electroweak penguin also contribute signif- 
icantly and may even cancel the gluon penguin contribution. There are large uncertainties 
from hadronic matrix evaluation, A^ CDS> dependce and errors in the CKM matrix. The 
range for e'/e is predicted to be between 10 -4 to 10~ 3 for A4 C = 300 MeV Hl4| . This is 



consisten with the experimental constraints from Fermilab, (7.4 ± 6.0) x 10 4 and CERN, 
(23 ±6.5) x 10- 4 pf. 

In the S3 x Z3 model there are several contributions to e'/ e. Due to large neutral Higgs 
masses, the neutral Higgs boson contributions to e'/ e are very small. However there may be 
large contributions from the charged Higgs bosons. The dominant contribution is from the 
charged Higgs boson mediated gluon penguin. The relevant AS* = 1 effective Larangian is 
given by 

L AS=1 = l da^Ch l Ap L + h l -^)^ a sG% , (19) 
where G a is the gluon field strength, A a are the SU(2>)c generators, and 

'h+ 

1 

I 



h = 32^M^ ( 2 " in ]^) Jro {(Vi f ^7«)iiWVicM7iy)B}C/ 

h + h 

3 

U = f^Xpr(| ~ In^ImiiVKMYf^Un'VKMlkMCf , (20) 



where Q = (a s (m/ l )/a s ( y u)) 14//23 ~ 0.17 is the QCD correction factor, and I is summed over 
u, c and t. We will use a s (fi) ~ 4n/6 for fi = lGeV. The above effective Lagrangian will 
generate a non-zero value for ImA fl6fl . Las=i also generates a non-zero value €ld for 
CP violation in K° — K° mixing due to long distance interactions through K° and tt, n, n' 



mixings [17]. One obtains [17,18 
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— w -0.196D , 

2m K ImM^ LD « 0.8 x l(r 7 (/i + / 2 )(GeV 3 ) , (21) 

where D is a supression factor of order 0(m 2 K , mfy/A 2 . ^/Ild is of order -0.014 to -0.1. 

We find that in both a) and b) cases, the donimant contributions are from the top quark 
in the loop arising from mixing in the charged Higgs boson couplings. For case a), we have 

e LD {ht) « o< , — i- fn-f , (22) 



m? + 150GeV m 2 



with a\ = 18, a2 = 25 and 03 = —7. 
And for case b), we have 



Uo » -7.35 x ^M^i)^^^ . (23) 



The contributions to e can be significant in both cases depending on the Higgs boson masses 
and the CP violating parameter 7771(722723) • We will study constraints on these parameters 
in Sec. VI. When these constraints are taken into account, €ld is generally constrained to be 
less than 3 x 10~ 5 for case a) and e'/e to be less than 3 x 10~ 5 . However, for case b), Eld 
can still be as large as 10~ 3 and e'/e can be 10~ 3 . 

V. CP VIOLATION IN THE NEUTRAL B SYSTEM. 

There are many processes which can test CP violation in the neutral B system. Some 
particularly interesting ones are Q 

B d -> J/ijKs , B d -> 7T+7T- , B s -> pK s . (24) 

The differences of time variation of decay rates for the above processes and their CP tran- 
formed states are given by 

_ r(B°(t) -> f C p) - T(B°(t) -» fcp) 

afCP r(B°(t) -> / p) + r(s°(t) - /op) 

(1 - \X\ 2 )cos(AM B t) - 2Im\sin(AM B t) , nr . 

= TTup ' (25) 
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where fcp indicates the final states. A is defined as 

where {q/p)B — \J Mf 2 * / M^, A and A are the decay amplitudes. If the final state contains 
K S) S = (q/p)K which has a phase of order 1CT 3 . For other cases S is equal to one. 

Non-zero asymmetry df C p signals CP violation. If |A| is not equal to one, it indicates 
that CP is violated in the decay amplitudes. In the MSM |A| is equal to one to a very good 
approximation for the above three processes. The asymmetries are proportional to ImX. In 
the MSM, the processes in Eq.(24) measure the three angles a, (3 and 7, 

ImX(B d -> J/ipKs) = -sin2f3 , 

Im\(B d -> 7r + 7r~) = sin2a , (27) 
ImX(B s — > pKs) = —sin2j , 

In the S3 x Z 3 model, the situation is very different. Although the CP violating decay 



amplitudes A and A are small, the phase of y M^*/M^ in the B — B mixing due to neutral 
Higgs boson exchange can be large. In case a), there is CP violation arising from the phase 
in Yukawa coupling of Higgs bosons, as well as CKM matrix, but the former is much larger. 
The three meaurements in Eq.(24) do not measure the angles a, f3 and 7 anymore. The first 
two processes will mostly measure the phases in M 12 d . We have 

Im\(B d -> 7T+7T-) ss Im\(B d -> J/^K S ) < 0.42 , from Rehl, 

ImX(B d -> 7T+7T-) ss JmA( J B d -> J/^K 5 ) < 0.19 , /rom i?e^ , (28) 

ImX{B d -> tt+tt") ss ImX{B d -> J/^K 5 ) w 0.19 , /orm i?e^ . 

For case b), we find 

imA(5 d -> 7T+7T-) w JmA( J B d -> J/^5) « -0.25 . (29) 

7mA for B s — > p-ft's is different for a) and b). For case a), the neutral Higgs boson 
contributions to the asymmetry are small. However ImX(B s — > pKs) due to CP violation 
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in the KM-matrix can be about 0.1. For case b), Im\(B s — > pKg) from neutral Higgs boson 
exchange is only about 0.02. 

If interpreted as in Eq.(27), we find for case a), sin2a = —sin2/3, sinj = 0.05, and 
a + (3 + 7 7^ 180°. For case b), we have, sin2a = —sin.2/3, siwy = 0.01. We again find, 
a + p + 1 ^ 180°. 



VI. THE NEUTRON AND ELECTRON ELECTRIC DIPOLE MOMENTS. 

The EDMs of neutron and electron in the MSM model are extremely small. The neutron 
EDM D n can only be generated at three loop level. It is predicted to be less than 10~ 31 ecm 
|19| . The electron EDM is even smaller (< 10~ 36 ecm) ||20|| . The experimental upper bound 
on the neutron EDM is 1.2 x 10 -25 ecm |2TH. For the electron the bound is about 10 -26 ecm 



22|| . If future measurement will obtain an EDM larger than the MSM model prediction, it 
will be an indication for new physics beyond the MSM. 

The prediction for the EDMs in the S3 x Z 3 are very different from the MSM. They may 
reach the experimental bounds. 

At the one loop level, the neutral Higgs contributions to the neutron EDM are small. 
For case a) we find that D n < 2 x 10 -28 ecm For case b), we have D n (d) ~ 2 x 10 _29 ecm. 
The u quark EDM is zero at the one loop level. 

There may be large contributions to the neutron EDM at the two loop level from the 



Weinberg operator [23| D n (W) and from the color dipole moment of gluon due to Bar-Zee 



type of diagrams f24l , |25l D n (BZ). In our model, we have 



1 22 
777 77? 

D n (W) w eC, w k-—ImZ\ t -l-ln ' ; 



64vr 2 "m 2 m 2 ' 

i i 
i \ i / 

where (w ~ 6 x 10~ 6 , and ( bz w 10~ 2 are the QCD correction factors, c u = 2 and q = 4, 
and A w lGeV is the chiral symmetry breaking scale. The parameters ImZ are given by 
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ImZ lt = ^2 Re (Yk,3 3 ( a ki ~ ia' ki ))Im(Y£ 33 (a kH - ia' kH )) , 

ImZ ^ = Z^ZT Im{ yk U A a ^ ~ i°hi)Yk>,u(<Xk'i + iu'k'i)) • (31) 

For case a), because there is no CP violation in the up quark sector only down quark loops 
contribute, D n (W) from the Weinberg operator at the two loop level is small. There are 
non-zero D n (BZ) from d-quark due to Bar-Zee mechanism. We find that the contributions 
from Refil l is also small (< 4 x 10~ 28 ecm). Reh® contribution is even smaller. 

For case b), the two loop contributions to the EDM are significantly larger because in 
this case there is CP violation in the top quark interaction. We have 

D n (BZ,u) « (2 ~ 8) x 10~ 26 ecm , 

D n (BZ, d) « (2 ~ 8) x 10- 27 ecm , (32) 

for m t between 100 GeV to 200 GeV. The contribution from the Weinberg operator is small, 
D n (W) < 10- 30 ecm. 

The charged Higgs bosons can also contribute to the neutron EDM. At the one loop 
level, the u and d quark EDM are given by 

d d = ^^ln^Im[^(V KM Yf) a (Vi^) u ] . (33) 

Z^i/T ill, _i_ ill, _i_ 

h- h- 

For d u , I is summed over d, s, and b; and for da, I is summed over u, c, and t. At the two 
loop level, there is a large contribution from the Weinberg operator, 

1 2 2 

I TYl Tfl 

D n (W) » eC^A— - 2 ImZH-f In-f , (34) 
6MX m h+ m h+ 



where ({y = 3 x 10 is the QCD correction factor, and 

ImZ* = —I m l lja l t (V KM Yf) 33 (V KM Y k u ) 33 ] . (35) 
mbm t 
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We find that in case a) the dominant contributions are from the two loop Weinberg 
operator. We have 



t GeV 2 , m 2 m 2 



D n (W) « h x lo-i«^-Zn ecm , (36) 

(150GW) 2 



where b x = 1.6, 6 2 = 1.4 and 6 3 = 1.2 x 10 -6 . 

Requiring the contributions to be less than the experimental value, we find the masses of 
have to be larger than 2.5TeV. There is no constraint on mass. Combining this in- 
formation with those from Eqs.(22) and (23), we find the charged Higgs boson contributions 
to Ild is less than 3 x 10~ 5 , and e'/e is less than 3 x 10 -5 . 

For case b), we find the dominant contribution is from the one loop d quark EDM. We 
have 

DM „ 5.4 x r-Mwi)^^ ecm . (37) 

Requiring D n (d) to be less than the experimental value, €ld is constrained to be less than 
10 -3 , and e'/e can still be of order 10~ 3 . Assuming maximum mixing, the mass of n + is 
constrained to be larger than 5 TeV. 

The S3 x Z 3 model may also have interesting CP violating signatures in the lepton 
sector. We assume the same S3 x Z3 assignments for the left handed and the charged right 
handed leptons as their quark partnenrs 0. The mass matrix and Yukawa couplings for 
the charged leptons are similar to the down quarks. One simply changes the parameters 
(a, b, c, d , and £) for quarks to (a; ,6; ,q , di , and £j = |£|e l<7 ') for leptons. We use [§: 
a\ = O.lOeGeV^ , b[ = , q = 1.781GeK , di = 8.6 x 10~ 3 GeV. For this set of parameters, we 
have m e = O.SllMeV^, m M = 106MeK and m T = 1784MeV which are in good agreement 
with experimental data. 

The calculation for the electron EDM is similar to the neutron EDM. For case a) we 
find that the one loop contributions are small (< 10 29 ecm) with a' = 80°. However the 
two loop contribution due to Bar-Zee mechanism p^j26[| can be as large as 1CT 27 ecm, for 
m t < 200GeV. For case b), we find that the one loop and two loop contributions are small 
(< 10 _33 ecm). 
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VII. CONCLUSIONS 



We have studied in detail some effects due to two different CP violating mechanisms 
in the S3 x Z 3 model. Both mechanisms discussed in this paper can explain the observed 
CP violation in the neutral K system. CP violation in the neutral K system and the mass 
difference in the neutral B system constrain the neutral Higgs boson masses to be in the 
multi TeV region. In the previous discussions we have chosen a particular set of parameters. 
The detailed predictions will depend on all the parameters, but the general features will 
remain to be the same. We have checked the predictions using another set of parameters, 
we find the changes are not significant except that the electron EDM for case b) can reach 
10~ 28 ecm. The predictions presented here represent the typical values for the observables. 
We summarize our results for e'/e, CP violation in B decays and the neutron and electron 
EDMs in Table 1. It is clear from Table 1, that the predictions for the observables considered 
are very different in the MSM and in multi-Higgs doublet models. Future experiments will 
be able to rule out some models. 



Observable 


MSM 


Case a) 


Case b) 


e 


Input 


Input 


Input 


e'/e 


io~ 4 ~ 10- 3 


~ 1(T 5 


~ 10~ 3 


B Decay 
Asymmetry 


a + + 7 = 180° 


a + (3 + 7 ^ 180° 
sin2y 7n —sin2f3 r* 0.2 ~ 0.4 
siny < 0.1 


a + (3 + 7 ^ 180° 
sin2y w —sin2f3 rs 0.25 
siwy ~ 


D n (ecm) 


1(T 31 ~ 10~ 33 


can reach 10~ 25 


can reach 10~ 25 


-D e (ecm) 


< 10~ 36 


< 10- 27 


< 10- 28 
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